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We present a rigorous analysis of the phenomenon of decoherence for general A^-level systems 
coupled to reservoirs of free massless bosonic fields. We apply our general results to the specific case 
of the qubit. Our approach does not involve master equation approximations and applies to a wide 
variety of systems which are not explicitly solvable. 
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I. INTRODUCTION 

We examine rigorously the phenomenon of quantum 
decoherence. This phenomenon is brought about by the 
interaction of a quantum system, called in what follows 
"the system S", with an environment, or "reservoir R". 
Decoherence is reflected in the temporal decay of off- 
diagonal elements of the reduced density matrix of the 
system in a given basis. The latter is determined by the 
measurement to be performed. To our knowledge, this 
phenomenon has been analyzed rigorously so far only for 
exphcitly solvable models, see e.g. [1-7]. In this paper we 
consider the decoherence phenomenon for quite general 
non-solvable models. Our analysis is based on the mod- 
ern theory of resonances for quantum statistical systems 
as developed in [8-15] (see also the book [16]), which is 
related to resonance theory in non-relativistic quantum 
electrodynamics [9, 17]. 

Let t) = [)s (8i f}R be the Hilbert space of the system 
interacting with the environment, and let 



be its Hamiltonian. Here, Hs and are the Hamilto- 
nians of the system and the reservoir, respectively, and 
Xv is an interaction with a coupling constant A € R. In 
the following we will omit trivial factors Is® and (8i11r. 
The reservoir is taken initially in an equilibrium state at 
some temperature T = 1//3 > 0. Let pt be the density 
matrix of the total system at time t. The reduced density 
matrix (of the system S) at time t is then formally given 

by 



Pt 



Tr 



(2) 
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where Tr^ is the partial trace with respect to the reser- 
voir degrees of freedom. Formulas (1) and (2) describe 
the situation where a state of the reservoir is given by a 
well-defined density matrix on the Hilbert space f)R. In 
order to describe decoherence and thermalization we need 
to consider "true" (dispersive) reservoirs, obtained for in- 
stance by taking a thermodynamic limit, or a continuous- 
mode limit. We refer to [18] for a detailed description of 
such reservoirs, which is not needed in the presentation 
of our results here. 

Let p(/3. A) be the equilibrium state of the interact- 
ing system at temperature T — 1/(3 and set p{f3,X) := 
TrRp(/?, A). There are three possible scenarios for the 
asymptotic behaviour of the reduced density matrix as 
t — > oo: 

(i) Pt — 'Poo =P(/3,A), 

(ii) Pt — ^ Poo 7^P(AA), 



(1) (iii) P( does not converge. 



The first situation is generic while the last two are not, 
although they are of interest, e.g. for energy conserving, 
or quantum non-demolition interactions characterized by 
[Hs,v] = 0, see [3, 18]. 

Decoherence is a basis-dependent notion. It is usu- 
ally defined as the vanishing of the off-diagonal elements 
[Pt]m, ri, m =/= n in the limit t — > oo, in a chosen ba- 
sis. Most often decoherence is defined w.r.t. the basis of 
eigenvectors of the system Hamiltonian Hg (the energy, 
or computational basis for a quantum register), though 
other bases, such as the position basis for a particle in a 
scattering medium [3], are also used. 

Since p(/9. A) is generically non-diagonal in the energy 
basis, the off-diagonal elements of will not vanish in 
the generic case, as t oo. Thus, strictly speaking, 
decoherence in this case should be defined as the decay 
(convergence) of the off-diagonals of to the correspond- 
ing off-diagonals of p(/3. A). The latter are 0(A). If these 
terms are neglected then decoherence manifests itself as a 
process in which initially coherent superpositions of basis 



2 



elements become incoherent statistical mixtures, 



III. QUBIT 



'^Cj,k\->Pj){'il}k\ — > '^Pj\'ipj){>Pj\, as t ^ oo. 

j,k j 

In particular, phase relations encoded in the Cj,fe disap- 
pear for large times. 



II. 



GENERAL RESULTS 



We consider A'^-dimensional quantum systems interact- 
ing with reservoirs of massless free quantum fields (pho- 
tons, phonons or other massless excitations) through an 
interaction v = G®ip{g), see also (1) and (6). Here, G is 
a hermitian N x N matrix and (p{g) is the bosonic field 
operator smoothed out with the form factor g{k), k gM.^. 
For any observable A of the system we set 



(^)t ■■= T^siPtA) = Trs+R(pt(^ ® Ir)). 



(3) 



Assuming certain regularity conditions on g{k) (allowing 
e.g. g{k) = \k\P e~^'^^"^ gi{a) where gi is a function on the 
sphere and where p= — 1/2 + n, n = 0, 1, . . ., m = 1, 2), 
we show that the ergodic averages 



((^»c 



lim 



{A),dt 



exist, i.e., that (^4)^ converges in the ergodic sense as 
i — > 00. Furthermore, we show that for f > 0, and for 
anyO<T'< ^, 

£5^0 

where the complex numbers e are the eigenvalues of a 
certain explicitly given operator K{t'), lying in the strip 
{z gC I < Im^; < r'/2}. They have the expansions 



£ = 4*) =e-A25W+0(A^), 



(5) 



where e G spec(i?s<8>]ls — lls<8>-ffs) = spec(ils) — spec(i?s) 

(s) 

and the Se are the eigenvalues of a matrix Ag, called a 
level-shift operator, acting on the eigenspace of Hg (8> Is — 
Is^Hs corresponding to the eigenvalue e (which is a sub- 
space of [)s $5 f}s)- The level shift operators play a central 
role in the ergodic theory of open quantum systems, see 
e.g. [18, 19]. 

The coefficients Re{A) in (4) are linear functionals 
of A which depend on the initial state po and the 
Hamiltonian H. They have the expansion Re{A) = 
J2{m,n)eie ^rn,nAm,n + 0(A2), where 4 is the collection 
of all pairs of indices such that e = Em — En, the Ek 
being the eigenvalues of Hg. Here, A^.n is the {m,n)- 
matrix element of the observable A in the energy basis 
of Hs, and the Km,n are coefficients depending on the 
initial state of the system (and on e, but not on A nor 
on A). 



Our results for the qubit can be summarized as follows. 
Consider a linear coupling, 



a c 
c b 



(S)(p{g), 



(6) 



where ip{g) is the Bose field operator as above. The form- 
factor g e I/^(M^,d^fc) contains an ultra-violet cutofi^ 
which introduces a time-scale tuv- This time scale de- 
pends on the physical system in question. We can think 
of it as coming from some frequency-cutoff determined 
by a characteristic length scale beyond which the inter- 
action decreases rapidly. For instance, for a phonon field 
Tuv is naturally identified with the inverse of the Debye 
frequency. We assume Tuv to be much smaller than the 
time scales considered here. 

A key role in the decoherence analysis is played by the 
infrared behaviour of form factors g{k). We characterize 
this behaviour by the unique p > — 1/2 satisfying 



< 



C < 00. 



(7) 



The power p depends on the physical model considered, 
e.g. for quantum-optical systems p = 1/2. We can treat 
p = -1/2 + n, n = 0, 1, . . .. 

Decoherence in models with interaction (6) with c = 
is considered in [1-6, 18, 21-23]. This is the situa- 
tion of a non-demolition (energy conserving) interaction, 
where [v, Hs] = 0, and consequently energy-exchange 
processes are suppressed. The resulting decoherence is 
called phase-decoherence. A particular model of phase- 
decoherence is obtained by the so-called position-position 
coupling, where the matrix in the interaction (6) is the 
Pauli matrix cr^ [2, 6, 22, 23]. On the other hand, energy- 
exchange processes, responsible for driving the system to 
equilibrium, have a probability proportional to jcj^", for 
some n > 1 (and a, b do not enter) [9, 10, 13, 15, 19, 20]. 
Thus the property c ^ is important for thermalization 
(return to equilibrium). 

We express the energy-exchange effectiveness by the 
function 



lim 



-/ 



d^/ccoth 



m 

2 



\9ik)\' 



■?7)2+e2' 



where rj > represents the energy at which processes 
between the qubit and the reservoir take place. Let 
A = E2 — E1 > be the energy gap of the qubit. In works 
on convergence to equilibrium it is usually assumed that 
|c|^^(A) > 0. This condition is called the "Fermi Golden 
Rule Condition". It means that the interaction induces 
second-order (A^) energy exchange processes at the Bohr 
frequency of the qubit (emission and absorption of reser- 
voir quanta). The condition c 7^ is actually necessary 
for thermalization while ^(A) > is not (higher order 
processes can drive the system to equilibrium). Observe 
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that £,{A) converges to a fixed function, as T — !■ 0, and 
increases exponentiaily as T ^ oo. The expression for 
decoherence times involves also ^(0), see (10). 

Our analysis allows to describe the dynamics of sys- 
tems which exhibit both thermalization and (phase) de- 
coherence. 

Let the initial density matrix, pt=Q, be of the form 
Po ® Pr,/3- (Our method does not require the initial state 
to be a product, see [18].) Denote by Pm.n the operator 
represented in the energy basis by the 2x2 matrix whose 
entries are zero, except the (n, m) entry which is one. We 
show that for t > 

[PtKi - ((f i,i»oo = [^0 + 0{X')] (8) 



and 



[Pt]l,2 - ((P2,l)), 



_ iteA(A) 



[Ca + 0{X')] (9) 



-l-e 



itso (A) 



0{X^) + e'*"-^(^)0(A2) + 0(A' e-*"'/2). 



Here, Co, Ca are explicit constants depending on the ini- 
tial condition Po J but not on A, and the resonance energies 
£ have the expansions 

£o(A) = iX\^\c\^^{A) + 0{X'') 

£a(A) = A + A^i^ + iAV [|c|2^(A) + (6-a)2e(0)] 

+0(A^) (10) 



and £_a(A) = — £a(A), with the real number 
R=^{b^-a^){g,w-'g) 



1 l^|2 
"2 



cl^P.V. / u^\g{\u\,a)\'^ coth 

MxS^ 



\ 2 J u-A' 



The error terms in (8), (9) and (10) satisfy (for small A): 



O(A^) 
A2 



< C and supj>o 



A^e- 



'/2 



< c. 



To our knowledge this is the first time that formulas 
for the decay of off-diagonal matrix elements of the re- 
duced density matrix are obtained for models which are 
not explicitly solvable, and without using uncontrolled 
master equation approximations (see e.g. [22] and refer- 
ences therein). 

Remarks. 1) The corresponding expressions for the 
matrix elements [/)j]2.2 and [pj]2.i art; obtained from the 
relations [pt]2,2 = 1 — [Pt]i,i (conservation of unit trace) 
and [pt]2,i = [pt]l,2 (hermiticity of pj. 

2) If the qubit is initially in one of the logic pure states 
Po = \'Pj){Vj\j where Hgifij = Ej^j, j = 1, 2, then we find 
Ca = 0, and Co = e'^^/^^g/SA _^ ^y:i/2 ^j. j ^ ^ ^^^^ 

Co = ef^^ie^^ + l)-3/2 for j = 2, sec [18]. 

3) To second order in A, the imaginary part of ea is 
increased by a term oc (6 — a)^^(O) only if p = —1/2, 
where p is defined in (7). For p> —1/2 we have ^(0) = 
and that contribution vanishes. For p < —1/2 we have 

m = oo. 



4) It is easy to sec that ^(A) and R contain purely 
quantum, vacuum fluctuation terms as well as thermal 
ones, while ^(0) is determined entirely by thermal fluc- 
tuations; it is proportional to /3~^ = T. 

5) The second order difference D, defined by 
Im£o(A) - ImeA(A) = X^D + 0{X'^), is D = 
iTT^ [|c|2^(A) - (6 - a)2^(0)] . For £> > the popula- 
tions converge to their limiting values faster than the off- 
diagonal matrix elements, as t ^ oo (coherence persists 
beyond thermalization of the populations). For D < 
the off-diagonal elements converge faster. If the inter- 
action matrix is diagonal (c = 0) then Z) < 0, if it is 
off-diagonal then D > 0. 

6) For energy-conserving interactions, c = 0, it follows 
that full decoherence occurs if and only if 6 7^ a and 
^(0) > 0. If either of these conditions are not satisfied 
then the off-diagonal matrix elements are purely oscilla- 
tory (while the populations are constant), sec also [18]. 

Illustration. Let the initial state of S be given by a 
coherent superposition in the energy basis, 



Po 



1 1 
1 1 



(11) 



We obtain the following expressions for the dynamics of 
the reduced matrix elements, for all t > 0: 



't\m,m 
[Pt\l,2 

[pth,i 



e-PEr. (-1)' 
1- — - 

+R 



■ tanh 



^it£o(A) 



m • 



1,2, 



1 pite-A(A) 



,iteA(A) 



+ RiA^,t), 
+ R2,i{X,t), 



where the numbers s are given in (10). The remainder 
terms satisfy \Rm,n{X,t)\ < CA^, uniformly in t > 0, 
and they can be decomposed into a sum of a constant 
part (in t) and a decaying one, Rm.n{X,t) = ((p„,,„))^ - 

^m,n^^+R'^,n{\ t), where li?;„,„(A, t)\ = 0(A2 e"^*), 
with 7 = min{Im£o, Im£±A}- Therefore, to second order 
in A, convergence of the populations to the equilibrium 
values (Gibbs law), and decoherence occur exponentially 
fast, with rates tt = [Im£o(A)]~-^ and tjj = [Im£A(A)]~-^, 
respectively. In particular, coherence of the initial state 
stays preserved on time scales of the order A~2[|c|2^(A)-|- 

{b-anm-\o.{. (10). 



IV. DISCUSSION 

Relation (4) gives a detailed picture of the dynamics 
of averages of observables. The resonance energies e and 
the functionals R^ can be calculated for concrete models, 
to arbitrary precision (in the sense of rigorous perturba- 
tion theory in A). See (8)-(10) for explicit expressions for 
the qubit, and the illustration above for an initially co- 
herent superposition given by (11). In the present work 
wc use relation (4) to discuss the processes of thermal- 
ization and decoherence of a qubit. In [18] we present, 
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besides a proof of (4), applications to energy- preserving 
(non-demolition) interactions and to registers of arbitrar- 
ily many qubits. It would be interesting to apply the 
techniques developed here to the analysis of the transi- 
tion from quantum to classical behaviour (see [1, 22]). 

In the absence of interaction (A = 0) wc have e = e € 
R, sec (5). Depending on the interaction each resonance 
energy e may migrate into the upper complex plane, or 
it may stay on the real axis, as A 7^ 0. The averages (^4)^ 
approach their ergodic means {{A))^ if and only if Ime > 

for all £ 7^ 0. In this case the convergence takes place 
on the time scale [Ime]~^. Otherwise {A)^ oscillates. A 

sufficient condition for decay is that Imd'f^ < (and A 
small, see (5)). 

There arc two kinds of processes which drive the de- 
cay: energy-exchange processes and energy preserving 
ones. The former are induced by interactions enabling 
processes of absorption and emission of field quanta with 
energies corresponding to the Bohr frequencies of S (this 
is the "Fermi Golden Rule Condition" [9, 13, 15, 19, 20]). 
Energy preserving interactions suppress such processes, 
allowing only for a phase change of the system during the 
evolution ("phase damping", [1-6, 21]). 

Even if the initial density matrix, pt=o, is a prod- 
uct of the system and reservoir density matrices, the 
density matrix, pt, at any subsequent moment of time 

1 > is not of the product form. The evolution cre- 
ates the system-reservoir entanglement. We develop a 
formula for {A)^ — {{A))^ for all observables A of any 
A^- level system S in [18]. If the system has the prop- 
erty of return to equilibrium, i.e., if ^(A) > 0, then 



in 



-/3E„ 



0(A2). 



Hence the Gibbs 
00 and then 



ooJm.n — "m,n Trs( c^'*"e 

distribution is obtained by first letting t 
A — » 0. A similar observation in the setting of the quan- 
tum Langevin equation has been made in [24]. If po is 
an arbitrary initial density matrix on Hs (S) Ti-K then our 
method yields a similar result, see [18]. 

Equations (8), (9) and (10) define the dccohcrcncc time 
scale, td = [Ini£A(A)]~^, and the thermalization time 
scale, tt = [Imeo(A)]~^. We should compare tjj with 
the decoherence time scales and with computational time 
scales in real systems. The former vary from 10*s for nu- 
clear spins in paramagnetic atoms to 10~^^s for electron- 
hole excitations in bulk semiconductors (see e.g. [25]). 

In the ubiquitous spin-boson model [26] , obtained as a 
two-state truncation of a double-well system or an atom 
interacting with a Bose field, the Hamiltonian is given 
by (1) with Hs = —^HAoax + iecr^ and v = ® fig)- 
Here, (7^, fz are Pauli spin matrices, e is the "bias" of the 
asymmetric double well, and Aq is the "bare tunneling 
matrix element". In the canonical basis, whose vectors 
represent the states of the system localized in the left and 
the right well, Hs has the representation 



Hs = 



1 



e -fiAo 
-hAo -e 



(12) 



The diagonalizati on of iJg yields Hg = diag(-E'+, i?_), 
where E± = ±^\/e^ + fi^Ag. The operator v = (Jz®^{g) 
is represented in the basis diagonalizing iJg as (6), with 



6 = -(^ + 1)-V2 and c = i(^ + 1)-' 



/2 
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